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摘 X. 随 着 植保 无 人 机 在 精细 农业 上 的 应 用 日 益 增 长 ， 目 前 在 植保 无 人 机 下 洗 风 场 演化 及 其 作用 下 的 雾 
滴 沉 积 膨 移 过 程 的 数值 模拟 方法 取得 了 快速 多 样 化 发 展 ， 但 对 各 方法 的 优势 、 缺 隐 、 适 用 范围 及 验证 手段 
仍 缺 乏 系统 的 梳理 。 本 文 针 对 无 粘 模型 、 计 算 流 体力 学 模型 及 格子 玻 尔 效 曼 模型 分 别 开 展 论述 。 基 于 涡 元 
法 的 无 粘 尾 调 模 型 优势 在 于 计算 过 程 简单 ， 但 由 于 缺乏 粘性 和 庙 流 模型 ， 其 雾 滴 沉积 款 移 模拟 精度 较 低 。 
计算 流体 力学 模型 又 分 为 有 限 体 积 法 与 有 限 差 分 法 。 其 中 ， 有 限 体积 法 鲁 棒 性 高 ， 可 适用 于 各 种 复杂 环境 
的 模拟 ， 但 格式 精度 有 限 ， 其 模拟 的 翼 尖 涡 耗 散 速度 远 快 于 实际 情况 ; 有 限 差 分 法 能 够 实现 对 辟 尖 涡 演 化 
的 高 时 空 精度 模拟 ,但 其 存在 网 格 结构 化 要 求 高 ， 算 力 要 求 过 大 等 问题 。 格 子 玻 尔 兹 曼 方 法 在 计算 具有 复 
杂 边界 条 件 和 非 平稳 运动 物体 的 三 维 流 场 问 题 中 具备 优势 ， 但 其 在 功能 多 样 性 和 完备 性 上 还 存在 不 足 。 上 
述 数 值 模型 精度 还 需 综 合 运 用 田间 实验 及 室内 实验 ， 如 高 速 粒子 图 像 测 速 (Particle Image Velocimetry, PIV) 
或 相位 多 普 勒 测速 (Phase Doppler Interferometry, PDI) 方法 进行 验证 和 优化 。 最 后 ， 本 文 提 出 了 未 来 植保 


无 人 机 施 药 模 拟 及 验证 方法 发 展 方向 。 
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verification methods of numerical modeling for plant protection unmanned aerial vehicle application[J]. Smart Agri- 
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=> 
1 引言 


不 同 于 地 面 作业 ， 农 业 航 空 植保 效果 受 环 境 
风 场 、 飞 机 飞行 高 度 和 速度 等 多 种 因素 影响 "”。 
根据 当前 环境 风 场 和 作业 参数 预测 雾 滴 球 移 沉积 


收 稿 日 期 : 2021-07-08 ”修订 日 期 2021-09-17 


区 域 ， 从 而 实时 调整 作业 参数 ， 是 农业 航空 精准 
化 植保 的 实现 手段 和 追求 目标 ”。 上 世纪 ?70 年 
代 ， 随 着 有 人 各 驶 固定 轰 飞 机 和 直升机 施 药 作业 
在 全 球 大 范围 普及 ,美国 国家 航空 航天 局 和 美国 
空军 开发 了 基于 拉 格 朗 日 方程 和 高 斯 模型 的 有 人 
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25h KELE 25 25 i8 SUE TR SEI CA Gricultural 
DISPersal, AGDISP), ， 目 前 已 成 为 全 球 通用 的 
雾 滴 球 移 预 测 工具 。 但 该 模型 缺乏 大 气 汕 流 直接 
模拟 能 力 ， 无 法 胜任 航空 施 药 幅 宽 内 复杂 气流 条 
件 作 用 下 雾 滴 运 动 沉积 的 准确 模拟 “。 

近年 来 ， 随 着 计算 机 技术 的 发 展 ， 直 接 基于 
计算 流体 力学 (Computational Fluid Dynamics, 
CFD) 数值 模拟 雾 滴 运 动 ， 尤 其 是 定 稼 流 场 状 态 
下 的 雾 滴 运动 规律 逐渐 成 为 可 能 。 目 前 研究 主要 
集中 于 地 面 风 送 式 喷雾 模拟 ， 风 场 、 喷 雾 结构 简 
JR. 未 涉及 包含 复杂 非 定 稼 旋 恤 风 场 作用 的 雾 滴 
运动 过 程 模拟 。 例 如 Tsay 等 * 对 无 冠 层 条 件 下 、 
Delele 等 ”和 Dekeyser 等 ”对 存在 冠 层 条 件 下 
的 空气 助力 喷雾 过 程 分 别 开 展 CFD 模拟 ， 评 佑 
了 不 同 作 业 参 数 下 的 空气 助力 喷雾 系统 效果 ， 从 
而 优化 作业 参数 、 评 佑 喷雾 球 移 特性 、 减 少 喷 雾 
DHA. Endalew 等 “" 对 空气 助力 地 面 喷雾 
机 开展 了 详细 的 数值 模拟 研究 ， 包 括 模型 构建 、 
气流 场 速 度 分 布 、 不 同 喷雾 机 型 影响 、 目 标 冠 层 
影响 等 。Duga 等 和 Baetens 等 基于 三 维 
CEFD 模 拟 了 风 送 雾 滴 的 球 移 运动 过 程 ， 考 虑 通过 
使 用 减少 飘移 和 标准 喷嘴 的 组 合 ， 避 人 免 地 面 沉 积 
损失 的 增加 ， 同 时 减少 药 液 飘 移 距 离 。 

随 着 植保 无 人 机 (Unmanned Aerial Vehicle, 
UAV) 在 精细 农业 上 的 应 用 日 益 增 长 ， 对 雾 滴 沉 
积 驹 移 过 程 的 模拟 面临 全 新 的 问题 。 以 当前 无 人 
机 植保 中 使 用 最 为 普遍 的 多 旋 己 飞机 及 无 人 直 升 
机 为 例 ， 在 超 低 空 飞行 植保 作业 时 ， 二 者 旋 忆 运 
动产 生 的 非 定常 下 洗 流 场 ， 具 有 周期 脉动 和 时 空 
分 布 不 均匀 的 特性 ， 与 固定 轰 有 人 驾驶 飞机 相对 
简单 的 二 维 化 尾 涡 结构 完全 不 同 ““ 。 二 者 施 药 
药 液 雾 滴 运动 受 低空 条 件 下 的 旋 波 痉 尖 涡 结构 演 
化 以 及 地 面 作 物 冠 层 及 旋 辟 下 洗 气 流 的 相互 作用 
影响 ， 预 测 其 沉积 飘移 规律 变 得 十 分 困难 。 

由 于 植保 无 人 机 施 药 应 用 呈现 爆发 式 增 长 ， 
对 施 药 效率 、 施 药 质量 、 环 境 污染 潜在 风险 评估 
的 需求 日 益 凸 显 ， 而 雾 滴 沉积 飘移 模拟 起 步 较 
晚 ， 面 临 研 究 难 度 大 ， 相 关 研 究 数 量 和 积累 又 远 


少 于 固定 路 有 人 机 施 药 研究 ， 无 法 及 时 满足 上 述 
需求 的 快速 增长 。 因 此 ， 本 文 拟 对 现 有 植保 无 人 
机 施 药 雾 滴 沉积 冉 移 数值 模拟 关键 技术 与 验证 方 
法 研究 进展 进行 梳理 ， 分 析 未 来 发 展 方向 ， 以 促 
进 该 领域 研究 取得 新 突破 。 

目前 关于 植保 无 人 机 施 药 研究 存在 的 主要 模 
型 研究 手段 包括 : 无 粘 模型 、 计 算 流体 力学 模型 
(包括 有 限 体积 法 和 有 限 差 分 法 模型 ) 以 及 格子 
玻 尔 兹 曼 模 型 。 本 文 就 这 些 模型 的 适用 范围 、 计 
算 开销 和 模拟 精度 等 方面 进行 了 总 结 ， 并 提出 了 
各 模型 未 来 优化 改进 研究 方向 。 本 文 还 针对 目前 
已 有 的 田间 和 室内 实验 数值 模型 验证 手段 提出 了 
发 展 优化 高 时 空 精度 的 非 定常 流体 测量 技术 ， 提 
升 模型 模拟 精度 的 建议 ， 为 植保 无 人 机 施 药 模 拟 
工作 发 展 提供 支撑 。 


2 植保 无 人 机 施 药 数值 建 模 关键 
技术 


2.1 无 粘 模型 


航空 施 药 雾 滴 肚 移 预 测 领域 最 为 常用 的 AG- 
DISP 软 件 即 建立 在 无 粘 方 程 基础 上 。 由 于 其 并 
不 直接 求解 纳 维 - 斯 托 克 期 (Navier-Stokes) 方 
程 ， 因 此 仅 能 针对 尾 涡 简 化 后 的 有 人 驾驶 直升机 
开展 雾 滴 球 移 模拟 ， 主 要 通过 将 固定 辟 或 有 人 轰 
驶 直升机 尾 涡 进 行 二 维 人 简化， 直接 模拟 二 维 涡 对 
近 地 发 展 过 程 。 虽 然 为 了 体现 粘性 影响 增加 了 尾 
涡 整 体 耗 散 系数 及 雾 滴 云 团 扩 散 系 数 ， 但 仍 无 法 
准确 模拟 植保 无 人 机 复杂 三 维 流 场 作用 下 的 雾 滴 
飘移 问题 。2018 年 ，Teske 等 "将 一 种 基于 涡 元 
法 的 无 粘 尾 涡 模 型 (Comprehensive Hierarchical 
Aeromechanics Rotorcraft Model, CHARM) 整合 
A dg HJ AGDISP 25 ij SUE THD EY, REDE T 
AGDISP Joys JF REGE Hb /] HY e ESE 76 ABL 25 25- 
T CERES TAI KY [eH o 
2.1.1 三 维 近 场 CHARM 模型 

基于 涡 元 法 的 无 粘 尾 涡 模 型 (CHARM) 于 
EHR 80 年 代 末 开发 完成 ， 其 主要 功能 是 模拟 
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尾 流 发 展 ， 最 初 用 于 人 研究 军用 直升机 旋 波 气流 对 
航母 甲板 人 员 的 有 影响， 军械 弹道 影响 ， 以 及 沙漠 
地 形 扬尘 影响 。 其 融合 并 改进 了 早期 算法 “ 旋 波 
飞行 器 (RotorCRAFT) ”和 “基于 影响 系数 的 
悬 停 表现 评估 (Evaluation of Hover Performance 
using Influence Coefficients, EHPIC) ”的 大 多 
数 功 能 "“”"， BLADE act ITE ELM. FEA 
非 升 力 面 的 源 / 偶 极 子 模拟 方法 ， 以 及 独 有 的 篆 
数 涡 量 云图 (Constant Vorticity Contour, CVC) 
全 恤 展 自由 尾 涡 模 型 mm 

(1) CVC BER. CVC 模型 能 够 用 于 模拟 复 
AR VER ER TAG EI. AN EEE AT 
速度 场 模 拟 ， 而 是 从 旋翼 旋转 平面 出 发 ， 基 于 流 
场 中 分 布 的 强度 不 变 的 大 量 涡 元 ， 将 其 按 当 地 涡 
量 分 布 进 行 跑 密 排 布 ， 并 依据 涡 元 诱导 速度 营 加 
来 获得 全 场 的 速度 分 布 。 涡 元 诱导 速度 遵循 毕 奥 
萨 伐 尔 定律 (Biot-Savart Law): 

"ES Dr,X ds 

2 4n ~ J 

EP, NARE; 为 位 置 拓 量 ; ds H 
ERE; 7” 为 到 涡 线 垂直 距离 ，m。CVC 模 型 的 
涡 线 自动 遵循 总 涡 量 守恒 的 刻 姆 堆 效 /开尔文 定 
律 ， 翼 尖 涡 旋转 半径 不 受 不 同 的 飞行 条 件 和 飞机 
型 号 影响 ， 但 存在 无 粘 点 涡 模 型 回 有 的 无 法 模拟 
尾 涡 耗 散 作用 和 庙 流 作用 的 问题 ”。 其 模拟 结 
果 如 图 1 所 示 。 


(1) 


ry 


图 1 基于 CHARM 模 拟 的 旋翼 尾 涡 演 化 过 程 5 
Fig. 1 Development of the rotor tip vortex simulated by 
CHARM simulation ^ 
(2) RRIA JI WIESO, eX 18767 E 77 II 
Bi ER HFE AE A Thee JE TR IA DLE f] 46 
Jb. HRT, NESEN h pu f e v SER ELA ib Tb 
组 成 的 四 边 形 表示 。 通 过 假设 每 个 旋 贺 四边形 控 


LS EAC PRE, ORF Pat Sc 
HERK. Kma PE a cs RE, RK 
fj EBA Mit. Ra DE A a BEE SE (Jou- 
kowsky Theorem) AE SE, R ff Jj R mi 7J i 
TAFE 

(3) 升力 / 非 升 力 面 的 源 / 偶 极 子 模型 。 对 于 
机 吴 表 面 引 起 的 涡 量 分 布 问 题 ， 主 要 采用 升力 / 
非 升 力 面 的 源 / 偶 极 子 模 型 进行 模拟 2 。 源 强 
度 设 置 为 环境 流 的 法 向 分 量 ， 然 后 调整 偶 极 子 强 
度 以 在 表面 产生 零 内 部 扰动 势 拉 。 通 过 对 表面 
势 值 进行 有 限 差 分 来 获得 表面 速度 和 压力 ， 通 过 
估计 分 离 点 位 置 来 确定 尾 流 平面 位 置 。 

(4) 分 级 涡 元 法 模型 。CHARM 开发 完成 后 
多 次 被 用 于 对 整 机 模型 悬 停 及 前 飞 状态 下 的 精准 
模拟 ， 并 逐渐 发 展 出 针对 多 旋翼 尾 流 模拟 的 功 
BEC". J T HEFt CHARM 模拟 的 实时 性 ， 研 究 
人 员 引 入 了 分 级 涡 元 法 (Hierarchical Fast Vor- 
tex，HFV)。 该 方法 最 初 基于 Appel ^ 开发 ,在 
计算 具有 NN 个 涡 元 的 自由 涡流 尾 流 问 题 时 ， 中 央 
Ab PEAS (Central Processing Unit, CPU) 每 次 迭 
代 计 算 量 增长 量 级 达 O (N’)。HFV 方 法 进一步 
将 计算 域 分 成 不 同 的 长 方 体 区 域 ， 并 基于 影响 系 
数 、 多 极 展开 和 泰勒 级 数 展开 ， 采 用 八 又 树 网 络 
将 不 同 区 域内 涡 元 进行 缩 并 ， 将 CPU 每 次 迭代 
计算 量 量 级 从 O N?) 降低 至 O (NlogN) ”。 
对 于 一 个 拥有 35,000 涡 元 的 算 例 而 言 ，HFV 方 法 
能 够 将 计算 量 减少 至 原始 算法 的 1%， 从 而 极 大 
提升 模拟 速度 。 
2.1.2 AGDISP Z ij SUE s 8H 

CHARM f5 93 3: 22 E] Ze ESE = AE i35 v6. AY 
的 速度 场 模拟 ， 在 远 场 速度 模拟 时 因 缺 乏 粘 性 和 
清流 模型 ， 模 拟 精 度 会 大 幅 下 降 ， 也 未 涉及 雾 滴 
运动 过 程 模 拟 。 因 此 作为 雾 滴 瞻 移 模型 ，AG- 
DISP 能够 与 其 实现 优势 互补 ， 共 同 为 植保 无 人 
机 施 药 模 拟 提供 完整 解决 方案 。AGDISP 采 用 的 
雾 滴 和 运动 模型 较为 简单 ， 主 要 包括 雾 滴 在 近 场 、 
远 场 和 冠 层 中 的 运动 模拟 。 

(1) 近 场 模拟 。 在 1000 m 区 域内 的 雾 滴 球 
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移 沉积 模拟 方面 ，AGDISP 采 用 基于 拉 格 朗 日 方 
法 的 雾 滴 运动 模型 ， 主 要 考虑 雾 滴 重 力 及 气动 忠 
力 对 雾 滴 运动 的 影响 ， 忽 略 雾 滴 附加 质量 力 、 
Saffman JJ, mMM EKRE, EJE K 


雾 滴 受 力 控制 方程 为 : 
d — t 
ga c-r) (2) 


HP, XNE, m; (U, - Vi) ATR 
EEZ, ms; 8 为 重力 加 速度 ，m/s; c, JR 
弛 豫 时 间 ，s。 

4p,D 
3Cyp,(U, - V.) 

其 中 ，p, 为 液体 密度 ，kg/m’'; p, 为 空气 密 
RE, kg/m; DAWA, m; Cj Cul 
阻力 系数 。 

上 述 控制 方程 最 大 程度 地 保留 了 雾 滴 空 间 中 
的 受 力 情况 ， 并 简化 了 雾 滴 和 运动 模拟 的 计算 量 。 
但 由 于 对 空气 速度 的 模拟 方面 缺乏 消 流 模型 的 文 
持 ， 其 无 法 对 大 气 边界 层 速度 脉动 量 和 注 动 能 


(3) 


Tp 


散 率 等 系列 参数 进行 准确 估计 ， 因 此 空气 速度 


的 估计 误差 会 被 引入 雾 滴 运动 受 力 方程 ， 进 而 影 
响 雾 滴 运动 模拟 准确 性 。 

(2) 远 场 模拟 。 为 了 减 小 计算 量 ，AGDISP 
在 1000 m 外 的 远 场 区 域 不 再 计算 具体 单个 雾 滴 
运动 过 程 ， 而 是 采用 基于 高 斯 分 布 的 雾 滴 云 团 分 
布 估计 ， 雾 滴 体 积 空间 分 布 为 : 


_ fe — BIEN 
V | SIE 


xtx,-bx)v : 
puo EUN jq HT 
_ u _ ü 
M, m exp| -5 " (5) 
x+x,- bx)v 
ioma - ) 
im Xoo 
Í 2jH,, =H ee 
- u (6) 
exp E ? 


(7) 


其 中 ，/ 为 雾 滴 体 积分 数 ; ONE. m/s; 
Xx 为 下 风 距 离 ，m; y 为 流 问 距离 ，m; z 为 高 度 ， 
m; 五 为 作业 高 度 ，m; 也 ,为 大 气 边界 层 混合 高 
E, m; v 为 雾 滴 沉降 速度 ，m/s; 五 为 平均 风速 ， 
m/s; 0,、0; 为 雾 滴 云 团 流向 、 垂 直 标 准 差 ; j 为 
用 于 加 和 的 系数 1 和 2。 上 述 方程 引入 了 大 气 边 
界 层 混 合 高 度 ， 将 大 气 稳定 性 的 部 分 影响 纳入 雾 
滴 沉 积 预 测 模型 ， 因 此 在 简化 计算 量 的 同时 具有 
良好 的 实用 性 。 

(3) 冠 层 模拟 。 冠 层 是 植保 无 人 机 施 药 作业 
模拟 的 重要 影响 因素 之 一 ， 由 于 冠 层 结 构 有 复杂 
性 和 非 均 匀 性 ， 其 模拟 难度 也 很 高 。AGDISP 中 
的 冠 层 模型 做 了 二 维 人 简化 ， 忽 略 冠 层 实际 结构 ， 
采用 不 同 作物 冠 层 类 型 的 穿 透 概率 、 液 滴 的 路 
径 长 度 和 叶 面 持 留 效率 来 确定 树冠 拦截 率 。 雾 
滴 在 冠 层 内 部 的 沉积 过 程 是 首先 接触 上 层 叶 面 ， 
之 后 它 可 能 会 从 表面 反弹 或 破碎 成 许多 较 小 的 
液 滴 "“。 为 解决 这 一 问题 ，Schou 等 ”1 引入 
Attané 等 3 发 明 的 模型 ， 用 来 描述 液 滴 在 叶 面 
的 铺展 行为 模式 。 假 设 铺展 雾 滴 为 一 扁平 圆柱 ， 
其 控制 方程 如 公式 8): 


Laa, tira xw. 
| “(I= 
E la fe 


1 21 
十 一 | 十 4 十 一 一 十 
cos 0) =| 40h 3 35 


-2/3 dr : 
1.410h ^r]|—| =0 (8) 
dt 
其 中 ，7 为 圆柱 半径 ，mj; CHR], s; 0 为 


接触 角 ，(°*); Obnesorge & Oh = —“—; wh 
poD, 

动力 粘度 系数 ，N'smm; pNP, kg/m’; DW 

液 滴 直 径 ，m。 对 液 滴 叶 面 碰撞 破碎 过 程 ， 引 入 
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临界 碰撞 参数 KK。 
K = OhRe'* (9) 


其 中 ，Re = ye, 可 基于 Forster 的 理论 对 


临界 碰撞 参数 玉 进 行 估 计 ， 从 而 判定 雾 滴 是 否 
破碎 Un, 
2.1.3 WARENI 

综合 来 看 ，AGDISP (Ey ERRE H 09 25 38 9A 
移 预测 工具 具有 快速 、 稳 定 的 模拟 功能 ， 而 无 粘 
模型 CHARM 也 极 大 减 小 了 复杂 风 场 作用 下 的 雾 
滴 运 动 预测 计算 量 。 二 者 整合 后 的 模型 理论 上 能 
够 提供 无 人 机 作业 条 件 下 旋 波 复杂 下 洗 流 场 随时 
间 变 化 的 全 过 程 模拟 ， 以 及 在 此 气流 场 条 件 下 的 
雾 滴 实时 运动 过 程 模拟 。 该 模型 的 优势 在 于 
计算 过 程 简 单 ， 模 拟 耗 时 少 ， 但 由 于 缺乏 粘性 作 
用 和 应 流 模 型 ， 其 雾 滴 运动 模拟 精度 受到 影响 。 
2018 年 Teske 等 介绍 了 基于 CHARM 和 AGDISP 
模型 开展 的 针对 DP-12 5038 A ICON / Jj 38 Wj 
种 不 同 机 型 无 人 机 施 药 在 不 同 侧 风 条 件 下 的 雾 滴 
运动 过 程 的 模拟 结果 ， 标 志 着 植保 无 人 机 施 药 过 
程 模拟 拥有 了 完整 的 无 粘 模型 解决 方案 Uu 

Teske 32 H] CHARM 模型 模拟 了 植保 无 人 机 
尾 流 发 展 过 程 ， 之 后 将 AGDISP 模 型 中 的 雾 滴 运 
动 控制 方程 引入 ， 对 尾 流 作用 下 的 雾 滴 潜 在 运动 
过 程 进行 了 模拟 。 通 过 分 析 不 同 作 业 速 度 条 件 下 
的 下 洗 气 流 强 度 和 雾 滴 地 面 沉积 率 ， 提 出 了 植保 
无 人 机 施 药 临 界 速度 概念 ， 发 现 当 植保 无 人 机 作 
业 速 度 超 出 某 个 临界 速度 后 ， 由 于 下 洗 气 流 强度 
减弱 ， 施 药 沉 积 效 果 将 严重 下 降 。 

但 目前 CHARM 和 AGDISP 整合 后 的 模型 仅 
有 Teske 等 进行 了 简单 介绍 ， 其 实际 应 用 效果 仍 
不 明确 ， 也 远 未 完善 。 从 Teske 等 发 表 文章 的 情 
况 看 ， 首 先 AGDISP 的 风 场 模拟 部 分 是 准 二 维 投 
影 模式 ， 暂 时 未 与 CHARM 模 型 的 三 维 风 场 进行 
整合 。 因 此 目前 耦合 模型 只 能 沿用 CHARM 模拟 
的 三 维 风 场 并 在 其 中 加 入 AGDISP 8922 334 UE BE 
型 ， 从 而 导致 模拟 雾 滴 驹 移 距 离 受 限 ， 且 无 法 模 
拟 环 境 参数 变化 带 来 的 影响 。 其 次 AGDISP 中 包 
含 的 冠 层 模型 无 法 应 用 于 目前 的 耦合 模型 。 最 


后 ， 耦 合 模型 对 植保 无 人 机 施 药 雾 滴 沉 积 的 联合 
模拟 结果 尚未 得 到 任何 实验 验证 。 因 此 以 AG- 
DISP 和 CHARM 为 代表 的 植保 无 人 机 施 药 无 粘 
耦合 模型 还 需 更 多 田间 实验 对 比 测试 及 更 多 详细 
的 模拟 信息 披露 以 证 明和 完善 其 实用 性 。 


2.2 计算 流体 力学 模型 


随 着 计算 机 技术 的 快速 发 展 ， 基 于 CFD 的 
数值 模拟 方法 应 用 越 来 越 广泛 。CFD 方 法 可 选用 
多 种 满 流 模型 ， 能 够 准确 模拟 流体 粘性 和 潮流 耗 
散 过 程 ， 因 此 在 对 复杂 风 场 的 模拟 方面 相对 AG- 
DISP 软件 具有 天 然 的 优势 。 目 前 已 有 的 研究 主 
要 集中 在 无 人 机 悬 停 及 前 飞 条 件 下 的 风 场 分 布 及 
近 场 筋 滴 分 布 ， 以 定常 流 平均 场 模拟 为 主 ， 还 包 
括 地 面 、 冠 层 等 因素 对 风 场 的 影响 研究 。 但 由 于 
旋翼 尾 流 的 复杂 特性 ， 对 计算 资源 的 消耗 很 大 ， 
已 有 研究 尚 无 法 对 雾 滴 球 移 全 过 程 开 展 动 态 
模拟 。 

CED 的 主要 原理 是 通过 求解 基本 的 流动 控制 
方程 包括 连续 性 方程 、 动 量 守恒 方程 得 到 气流 速 
度 时 空 分 布 ， 其 控制 方程 如 下 。 

连续 性 方程 : 

V(pv) = 0 (10) 

Navier-Stokes 动 量 守恒 方程 : 


E *V(ov»xv)--Vp- Vi*pg (11) 

其 中 ，p 为 空气 密度 kg/m’; 7 为 气流 平均 
HERE, ms; p UR, N/m’; ?为 粘性 应 力 
张 量 ，N/m”。 由 于 潮流 的 复杂 性 ,目前 工程 上 通 
常 使 用 雷诺 平均 方程 (Reynolds-averaged Navier- 
Stokes equations, RANS) 处 理 Navier-Stokes 方 
程 。 雷 诺 平 均 将 漠 流 瞬时 运动 分 解 为 平均 运动 和 
脉动 运动 两 部 分 ， 把 脉动 运动 部 分 对 平均 运动 的 
贡献 通过 雷诺 应 力 项 来 模 化 ， 通 过 满 流 模型 来 封 
闭 雷 诺 平均 Navier-Stokes 方 程 。 除 了 流动 控制 方 
程 不 同 ， 计 算 流 体力 学 模型 中 的 雾 滴 运 动 控 制 方 
程 也 比 无 粘 方 法 要 复杂 得 多 ， 主 要 通过 对 拉 格 野 
日 坐标 系 下 颗粒 作用 力 的 微分 方程 贝 赛 特 - 鲍 瑟 
内 斯 克 - 奥 森 (Basset-Boussinesq-Ossen, BBO) 
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方程 进行 积分 来 求解 。 在 笛 卡 尔 坐标 系 下 ， 颗 粒 
的 BBO 方 程 为 : 


du, oP 
my 一 zi em ag 3nud, f (u; = u,) 十 
DAC = ù) " 
2 
3 Fc 
5d, nup: | ——— ds + 12 
54; nup, f. tm (m 


oP : pv, ú,- ú 
其 中 ， V, oy 是 压力 梯度 项 ; Ca 


ded 
vVt-s 
set 力 项 ; 3nud, f (u, - u,) Je ES CRESCE P] BEL 
" C;d,|u, - ,| "M 
力 项 ; 阻力 修正 系数 /一 uj 为 连 
续 相 的 速度 ，m/s; wu, 为 离散 相 上 颗粒 的 速度 ， 
m/s; /为 流体 的 动力 粘度 ，kg/ (ms); 2 为 流体 
WREE, kg/m; p, NBEO AE, kg/m; 
d, 为 颗粒 的 粒 径 ，m。 

液 滴 与 汗 流 的 相互 作用 则 通过 离散 随机 游 走 
模型 模拟 。 将 每 个 涡 定 义 为 高 斯 分 布 的 随机 速度 
波动 ， 即 ww'、v' 和 w'， 以 及 时 间 尺 度 t,。 速 度 脉 
ae MWA w=CV2k3, W ee HE A tt 
t, = 0.3k/e, FEA CA IEA A TH BALE, kA im ii 
动能 ，g 是 潮流 耗 散 率 。 

上 述 方程 一 般 通 过 离散 化 进行 数值 计算 求 
解 ， 具 体 方式 包括 有 限 体 积 法 和 有 限 差 分 法 等 。 
2.2.1 有 限 体 积 法 

由 于 上 述 流动 控制 方程 为 非 线 性 偏 微 分 方 
程 ， 无 法 得 到 精确 解 。CFD 通 过 数值 求解 得 到 这 
些 偏 微分 方程 的 近似 解 。 目 前 常用 的 计算 流体 力 
学 软件 Fluent、CFX、Openfoam 等 均 基 于 有 限 体 
积 法 离散 控制 方程 ， 将 偏 微分 方程 转换 为 代数 方 
程 组 ， 求 得 该 方程 组 的 解 ， 作 为 偏 微分 方程 的 数 
值 近似 解 。 

由 于 有 限 体 积 法 在 商用 计算 软件 中 的 普及 ， 
目前 针对 植保 无 人 机 的 风 场 模拟 结果 大 部 分 来 自 
有 限 体积 法 。 但 有 限 体 积 法 受 网 格 质 量 影响 较 
大 ， 其 常用 的 迎风 格式 精度 有 限 (2 阶 精度 )， 在 


—— : 
是 附加 质量 力 项 ; 7 d? map, 人 


ds 是 Bas- 


HI] bep m BERT P. EMERE SR iis AE AL 
限 差分 法 更 容易 产生 人 工 耗 散 ， 从 而 导致 其 模拟 
的 演 尖 涡 耗 散 速度 远 快 于 实际 情况 。 实 际 应 用 过 
程 中 容易 导致 因 计 算 资源 不 足 而 缩减 时 空 精度 或 
尺度 。 如 图 2 所 示 ， 低 精度 格式 产生 的 泪 尖 涡 在 
2 个 周期 内 耗 散 。 目 前 基于 商用 软件 进行 的 植保 
无 人 机 风 场 模拟 结果 表明 ， 有 限 体积 法 的 植保 无 
人 机 施 药 模 拟 尚 无 法 完整 复 现 旋 波光 尖 涡 发 展 过 
程 ， 其 非 定常 流 场 模拟 准确 性 受到 一 定制 约 。 


图 2 低 精度 格式 有 限 体积 法 模拟 时 楼 尖 涡 耗 散 过 快 
Fig. 2 Rotor tip vortex decay too fast when modeled by low 


order accurate FVM 


目前 关于 植保 无 人 机 下 洗 风 场 和 施 药 过 程 模 
拟 的 主流 研究 方向 为 有 限 体 积 法 ， 如 杨 风 波 
等 P" 基于 Fluent 模 拟 多 旋 慢 飞机 悬 停 状 态 下 的 
下 洗 气 流 场 速度 分 布 及 雾 滴 空 间 分 布 情况 ， 并 结 
合 室内 悬 停 实 验 进行 验证 。 模 拟 使 用 了 约 560 万 
个 网 格 ， 忽 略 了 机 身影 响 , 使 用 k-s 潮流 模型 。 
从 风 场 模拟 结果 看 ， 其 能 够 反映 整体 风 场 空间 分 
布 形 态 ， 标 记 点 平均 风速 模拟 值 与 实验 测量 结 
误差 在 9% 以 内 ， 具 有 良好 的 模拟 准确 性 ， 但 对 
瞬时 辟 尖 涡 精 细 结 构 发 展 的 模拟 损失 较 大 ， 难 以 
观察 到 翼 尖 涡 发 展演 化 过 程 。Shi 等 ^U 进行 了 前 
飞 状态 直升机 下 洗 风 场 和 雾 滴 喷 施 的 模拟 ， 使 用 
478 万 个 网 格 ， 采用 了 SST-k-o 消 流 模型 ， 成 功 
模拟 了 植保 无 人 机 施 药 筋 滴 地 面 沉 积分 布 。 该 和 雾 
滴 沉 积 模拟 结果 与 田间 实验 结果 吻合 度 良 好 ; 但 
对 风 场 发 展 并 未 关注 ， 难 以 判断 是 否 具有 良好 的 
风 场 结构 模拟 准确 度 。Zhu 等 ^U zr DLE 
下 洗 风 场 作用 下 的 雾 滴 沉 积 球 移 过 程 开 展 了 模 
拟 ， 使 用 22.8 万 个 网 格 ， 采 用 定常 模拟 旋 恤 下 洗 
风 场 。 之 后 将 粒子 注入 风 场 并 模拟 不 同 侧 风 条 件 
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下 的 雾 滴 偏 移 沉积 情况 。 该 模拟 未 关注 旋 驾 下 洗 
气流 的 非 定常 特性 及 波 尖 涡 演化 等 问题 。 张 宋 超 
等 ”模拟 了 N-3 型 农用 无 人 直升机 航空 施 药 雾 
滴 结 移 过 程 ， 采 用 130 万 个 网 格 以 及 面 源 喷 雾 模 


HA, FERRERA THER, EREE S 
期 后 仍 保持 基本 形态 。Lakshminarayan $ ^" jz 
用 高 精 格式 对 微型 共 轴 双 旋翼 悬 停 风 场 进 行 了 细 
致 模拟 ， 网 格 数 量 达 660 万 个 ， 重 点 关注 了 旋 要 


拟 ， 模 拟 结果 与 实验 结果 相差 较 大 ， 原 因 可 能 是 
实验 与 模拟 条 件 差 异 、 交 光 示 踪 剂 日 光 分 解 影响 
实验 结果 等 ， 模 拟 未 关注 风 场 分 布 情况 。 张 察 
Ak UU. MAE. CU MP EX F Ub CULTE IDE OR 
响 ， 以 及 冠 层 对 旋 豆 下 洗 气 流 的 影响 进行 了 模 
拟 ， 对 有 /无 冠 层 情况 分 别 采 用 734 万 和 3590 万 个 
网 格 进行 模拟 ， 与 实验 结果 对 比 最 大 模拟 误差 在 
20% EA, MERU KAZE A, HERRERA aa 
并 不 清晰 。 上 述 人 研究 主要 集中 于 对 植保 无 人 机 下 
洗 风 场 及 雾 滴 近 场 范围 沉积 的 模拟 ， 并 与 实验 结 
果 进 行 了 对 比 。 研 究 涉及 模型 结构 地 表 、 环 境 条 
件 多 样 ， 显 示 了 有 限 体积 法 良好 的 鲁 棒 性 和 适应 
性 。 但 受 有 限 体积 法 格式 精度 条 件 限 制 ， 对 流 场 
精细 结构 的 分 析 还 较 少 ， 也 尚未 涉及 大 范围 雾 滴 
运动 分 析 。 
2.2.2 有限 差分 法 

与 有 限 体积 法 不 同 ， 有 限 差 分 法 将 求解 区 域 
离散 成 差分 网 格 ， 以 有 限 个 节点 代替 连续 求解 
域 ， 用 网 格 上 的 差 商 代替 格 点 上 的 空间 导数 ， 从 
而 将 微分 方程 转化 为 以 网 格 节 点 参数 作为 未 知 量 
的 差分 方程 组 ， 然 后 进行 求解 。 

相对 于 有 限 体 积 法 来 说 ， 有 限 差分 法 所 采用 
的 结构 化 网 格 更 易于 构造 高 阶 精 度 的 数值 格式 。 
高 阶 数值 格式 如 加 权 本 质 无 震荡 (Weighted Es- 
sentially Non-Oscillatory, WENO) /本 质 无 震荡 
(Essentially Non-Oscillatory，ENO) 格式 等 可 达 
4~5 阶 精度 ， 配 合 自 适应 网 格 技术 ， 达 到 对 疲 尖 
涡 演化 的 高 时 空 精度 模拟 ， 能 够 复 现 真 实 旋 浑 下 
洗 流 场 的 典型 流动 结构 发 展 过 程 。 

目前 有 限 差 分 法 还 主要 应 用 于 基础 研究 领 
W, Ul Xu Fl Weng "" 基于 NACA0012 3 HY è Ey 
旋翼 ， 构 建 基于 不 同 高 精 格式 的 验证 模型 ， 使 用 
137 万 个 网 格 。 模 拟 结果 可 以 明显 看 出 高 精 格式 
在 相同 网 格 条 件 下 能 够 保留 更 多 的 旋 波 波 尖 涡 精 


访 面 流动 参数 模拟 精度 及 波 尖 涡 与 旋 波 作用 及 演 
化 过 程 的 模拟 准确 度 。Kalra ” Sl Fl RS AS 
J Be Te BE BER Ts Sa, LA eG WENO 格式 计算 
Bese PEK, OLA A BCI 2160 万 个 。 
研究 发 现 高 精 格式 算法 可 模拟 多 个 旋转 周期 后 的 
疲 尖 涡 结构 演化 过 程 ， 流 动 结构 与 室内 实验 结 
几乎 完全 吻合 ， 近 地 面 展 向 速度 剖面 仍 与 实验 结 
果 存 在 约 25% 的 误差 。 上 述 研究 相对 有 限 体积 法 
拥有 更 高 的 模拟 精度 ， 能 够 更 好 地 模拟 旋 注 吕 尖 
涡 运 动 和 耗 散 过 程 ， 但 都 未 模拟 复杂 机 呈 结 构 外 
形 ， 也 未 涉及 作物 冠 层 及 雾 滴 和 运动 过 程 模拟 ， 无 
法 为 植保 无 人 机 施 药 过 程 模拟 提供 直接 解决 
方案 。 

目前 有 限 差 分 法 在 实际 应 用 中 仍 面临 许多 问 
题 ， 如 对 于 复杂 外 形 的 模型 难以 构建 与 该 算法 相 
适应 的 结构 化 网 格 ， 算 法 存在 对 于 两 相 流 算法 的 
兼容 性 问题 ， 对 作物 冠 层 的 建 模 方 法 尚未 完善 ， 
对 大 空间 三 维 建 模 的 算 力 要 求 过 大 等 ， 目 前 仅 适 
用 于 简单 外 形 及 环境 参数 作用 下 的 流动 机 理 研 
究 。 基 于 有 限 差分 法 的 无 人 机 施 药 数值 模型 短期 
内 仍 难以 推广 。 


2.3 格子 玻 尔 兹 曼 模型 


近年 来 ， 一 种 基于 格子 玻 尔 效 曼 方法 (Lat- 
tice-Boltzman-Method, LBM) 的 计算 流体 动力 
学 方法 在 计算 具有 复杂 边界 条 件 和 非 平 稳 运 动物 
体 的 三 维 流 场 问题 中 逐渐 显示 优势 “1。 
2.3.1 LBM 模型 原理 

LBM 是 一 种 笛 卡 尔 坐 标 系 下 离散 点 对 应 速 
度 在 格子 节点 空间 分 布 的 粒子 化 方法 。 与 计算 流 
体力 学 模型 基于 连续 介质 假设 下 的 Navier-Stokes 
方程 不 同 ，LBM 基于 分 子 运动 理论 的 玻 尔 效 曼 
方程 ， 在 宏观 层面 兼容 Navier-Stokes 方 程 ， 并 拥 
有 更 广泛 的 适用 性 ， 主 要 通过 流 场 平衡 态 分布 函 
数 和 碰撞 算 子 进行 迭代 计算 。 一 种 典型 的 3 维 27 
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自由 度 (D3Q27) 的 LBM 如 图 3 和 所 示 。 


图 3 D3Q27 离散 速度 格子 模型 
Fig. 3 D3Q27 lattice model for velocity discretization 


此 方法 基于 粒子 碰撞 原理 ， 每 个 离散 单元 具 
备 较 高 自由 度 ， 整 体 自然 达到 4 阶 空间 离散 精 
度 ”， 能 够 以 高 精度 模拟 旋翼 辟 尖 涡 产 生 及 发 
展 过 程 ， 在 相同 网 格 数 的 前 提 下 比 传统 有 限 体积 
方法 拥有 明显 模拟 精度 优势 。 同 时 该 方法 的 并 行 
计算 效率 几乎 不 受 核心 数量 限制 ， 更 适合 开展 并 
行 计算 ， 提 高 模拟 速度 。 图 3 所 示 的 D3Q27 离 散 
速度 模型 参数 见 表 1， 其 中 声速 为 cv， 离 散 速度 
efa=0,1,...， 261, BUE w, c= 6x/5t， 其 中 
Ox K ot NAB FIN ZS AP KK 


X1 D3Q27 离散 速度 模型 参数 
Table 1 Parameters of the D3Q27 discrete velocity model 


模型 格子 声速 cs/c 格子 离散 速度 e,/c Ew, 
(0, 0, 0) 8/27 (a = 0) 
(41, 0,0), (0, +1, 0), (0, 0, +1) 2/27 (a = 1, +++, 6) 
D3Q27 /V3 
(x1, +1, 0), (x1, 0, £1), (0, x1, +1) 1/54 (a = 7, +++, 18) 
(41, £1, +1) 1/216 (a = 19, +++, 26) 


LBM 3E 3f iE 2-46 AS PREY AN E TG ET 
来 模拟 流 场 各 参数 的 演变 。 流 体 平衡 态 分 布 函 数 
f 的 发 展 方程 为 : 

f(r, + e,ôt,t + ót) - f (r, t) = 

-M^$[m(r,.:) - m" (r,t)|+ F (r,t) (13) 

其 中 ，M 是 动量 分 布 函数 的 线性 转化 矩阵 : 
m- Mf; 广 是 位 置 矢量 ; SEEE; mH 
动量 ，kg'm/s; FAS, kg:m/s; 表示 如 下 
矢量 : 

F (rat) = (Flen) f ret) (rr) (14) 

等 效 动量 m? = Mf, MARERE S 9927 X 
27 矩阵 。 
2.3.2 ”LBM 模型 特性 

LBM 模 型 自身 具备 一 些 独 特 的 特性 ， 有 利 
于 其 处 理 复杂 外 形 、 复 杂 流 场 条 件 下 的 高 精度 模 
拟 问 题 。 由 于 格子 模型 结构 符合 八 又 树 结构 ， 天 
然 具 备 自 适应 的 优势 ， 有 利于 处 理 复杂 外 形 、 复 
杂 流 场 等 的 计算 ， 既 可 以 对 旋翼 翼 面 、 机 身 、 冠 
层 等 复杂 表面 进行 加 密 ， 提 升 壁面 边界 层 模 拟 准 
确 性 ， 也 可 以 基于 涡 量 分 布 对 涡 量 集中 区 域 的 网 


格 尺 寸 进行 实时 加 密 ， 从 而 动态 提高 旋 恤 尾 涡 结 
构 状 识 精 度 。 因 此 更 有 利于 针对 植保 无 人 机 旋 滤 
下 洗 风 场 进行 高 精度 模拟 。 

基于 LBM 的 商用 软件 XFlow 采 用 壁面 适应 
局 地 涡 粘 模 型 。 由 于 LBM 受 格子 形状 所 限 ， 一 
般 壁 面 附近 的 网 格 尺寸 较 大 ， 需 要 匹配 专门 的 庙 
流 模型 以 提升 壁面 边界 层 计算 精度 ， 降 低 计算 
量 。 由 于 速度 梯度 张 量 在 当地 自动 满足 LBM， 
从 而 极 大 提升 了 壁面 大 涡 模拟 的 效率 。 而 对 有 限 
体积 法 来 说 ， 速 度 梯度 张 量 作为 空间 导数 ， 需 要 
临近 节点 的 信息 参与 计算 。 因 此 大 涡 模拟 壁 湛 流 
模型 在 LBM 中 相对 有 限 体积 法 拥有 更 高 的 运算 
效率 。 

针对 中 心动 量 空间 散射 模型 采用 多 松弛 时 间 
方法 ， 可 独立 调节 不 同 松弛 时 间 来 控制 稳定 性 ， 
相对 传统 的 Bhatnagar-Gross-Krook (BGK) 模型 
来 说 具有 更 高 的 稳定 度 。 另 外 还 克服 了 一 些 
BGK 模型 存在 的 明显 问题 ， 例 如 Prandtl 数 以 及 
动力 粘度 与 第 二 粘度 比值 无 法 改变 。 
2.8.8 LBM 模型 应 用 

目前 植保 无 人 机 施 药 研究 领域 ，LBM 已 获 
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得 初步 应 用 。 例 如 Zhang 等 中 运用 LBM 模 拟 了 
旋翼 下 洗 风 场 在 不 同 高 度 及 作业 速度 条 件 下 的 形 
态 变 化 ， 获 取 了 流 场 精细 结构 的 空间 分 布 演化 情 
况 。 与 田间 风 场 测量 结果 对 照 显示 ， 模 型 成 功 复 
Rr be ARE 显示 了 该 模型 良 
好 的 流 场 结构 模拟 准确 度 。 文 晟 等 “” 运用 
LBM 对 植保 无 人 直升机 及 多 旋 恤 植保 无 人 机 下 
洗 风 场 分 别 开 展 模拟 ， 其 翼 尖 涡 可 保持 约 三 个 周 
期 ， 还 显示 了 尾 桨 涡 和 螺旋 尾 涡 的 发 展 过 程 。 同 
时 在 风 场 中 按 喷头 实测 结果 加 入 不 同 粒 径 分 布 的 
雾 滴 ， 模 拟 了 雾 滴 在 不 同 作业 参数 下 的 近 场 分 布 
情况 。 但 上 述 研究 主要 集中 在 机 体 尾 流风 场 及 
幅 宽 内 却 滴 空间 分 布 上 ， 但 仍 存在 缺乏 与 相同 


条 件 下 的 实测 结果 详细 对 照 ， 缺 乏 远 场 雾 滴 运 动 
模拟 数据 等 问题 。 作 者 团队 “基于 无 人 机 实 
际 作业 参数 进行 全 尺度 三 维 空间 建 模 ， 以 地 面 
为 固定 参考 系 ， 实 现 了 无 人 机 巡航 状态 下 洗 风 
场 动态 建 模 。 模 拟 时 长 达 12 s 以 上 ,空间 范围 达 
80m (K) X60m ( 宽 ) X10m (高 )。 

图 4 显示 了 基于 LBM 模拟 获得 的 植保 无 人 
机 旋 恤 下 洗 速度 在 不 同 展 向 距离 上 的 剖面 与 田间 
实验 和 室内 粒子 图 像 测速 
结果 的 对 比 。 从 图 4 可 以 看 出 ， 
LBM 模 拟 结果 在 近 地 面 范围 相对 实验 结果 差距 
较 大 ， 在 远离 地 面 范 围 内 速度 模拟 精度 可 以 
接受 。 


(Particle Image Velo- 


cimetry, PIV) 
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图 4 植保 无 人 机 号 和 停 状态 下 各 展 向 蕉 面 速度 训 面 模拟 结果 与 田间 实验 及 室内 PIV 实验 对 比 


Fig. 4 Comparison of the simulation results of the speed of each spanwise section profile of the plant protection unmanned 


aerial vehicle in the hovering state with the field test and the indoor PIV test 


在 获得 流 场 运动 规律 的 基础 上 ， 对 拉 格 朗 日 
坐标 系 下 颗粒 作用 力 的 微分 方程 进行 积分 来 求解 
离散 相 颗 粒 的 运动 轨迹 ， 可 实现 雾 滴 沉积 规律 的 
预测 。 目 前 基于 LBM 的 商用 软件 XFlow 对 于 雾 
滴 运 动 建 模 的 功能 还 相对 简单 ， 雾 滴 运 动 模 型 采 
用 了 类 似 AGDISP 所 使 用 的 雾 滴 运动 方程 : 

dU, 

ro 50-5 tg, (15) 

其 中 ，U 为 雾 滴 速度 ，m/s; ao NE 
滴 气 动 阻力 造成 的 加 速度 ，m/s’; gi; 为 重力 加 
WE, m’. HTAR (15) REBAR 
型 ， 而 无 人 机 施 药 过 程 中 多 采用 高 浓度 药 液 ， 因 
雾 滴 蒸发 导致 的 沉积 膨 移 模拟 误差 相对 有 人 机 作 
业 模 拟 误差 较 小 ， 因 此 该 方法 目前 更 适合 用 于 无 


人 机 施 药 过 程 的 模拟 。 

一 组 针对 不 同 作业 高 度 的 植保 无 人 机 施 药 空 
间 分 布 情况 如 图 $ 所 示 。 从 图 中 可 以 看 出 ， 运 用 
LBM fj fe s H ES 7L d 3s y B) 
时 间 累 积 过 程 ， 获 取 不 同 作 业 状 态 下 不 同时 刻 的 
25 ii ax UAR 并 具备 较 高 的 时 空 分 辨 率 。 
这 一 优势 有 助 于 对 不 同 作 业 机 型 、 环 境 参数 、 喷 
雾 条 件 下 的 雾 滴 沉 积 球 移 特性 进行 快速 模拟 并 总 
结 出 相关 定量 性 规律 ， 从 而 为 分 析 优 化 作业 参数 
KALE HE EA HI SCF. FET LBM 模拟 植保 
无 人 机 施 药 雾 滴 地 面 沉 积分 布 结果 与 田间 实验 结 
果 对 比如 图 6 所 示 。 其 结果 显示 幅 宽 内 雾 滴 相对 
沉积 量 的 模型 预测 结果 和 田间 实验 结果 差异 在 可 
接受 范围 ，LBM 雾 滴 沉 积 模拟 结果 可 信 。 


10 
Ga oom 50m 40m 30m 20m lii t TERIS "J i 
(a) 作 业 高 度 为 1.5m 时 雾 滴 空间 分 布 (b) 作 业 高 度 为 1.5 m 时 和 雾 滴 运 动 范 围 边界 拟 合 
y=-0.17*x? -0.2*x + 60 E UN fing 
-20 -10 mo xx 10 20 30 
〈《c) 作 业 高 度 为 2.53m 时 雾 滴 空间 分 布 (d) 作 业 高 度 为 2.5 m 时 和 雾 滴 运 动 范 围 边界 拟 合 
=20 to 6 10 20 30 
展 向 距离 x /m 
(e) 作 业 高 度 为 3.5 m 时 和 雾 滴 空间 分 布 (人 作业 高 度 为 3.5 m 时 雾 滴 运 动 范围 边界 拟 合 
图 5 基于 格子 玻 尔 褒 曼 法 模拟 的 不 同 作业 高 度 植保 无 人 机 施 药 空间 分 布 规律 9 
Fig. 5 Droplet spatial distributions of UAV spraying for different flight heights based on LBM simulation 
米 A * 
a 3 数值 模型 的 验证 
90F W Numerical simulation by Tang et al. (2020b)| ^ 
gil: -8-Field test result by Chen et al. (2017) | 
eR x S ee ` y NARI 
» 上 述 无 粘 模型 、 计 算 流体 力学 模型 及 格子 玻 
E 尔 兹 曼 模型 模拟 结果 的 可 靠 性 和 准确 性 均 需 要 通 
& 过 对 植保 无 人 机 旋翼 下 洗 风 场 ， 施 药 雾 滴 地 面 沉 


图 6 基于 格子 玻 尔 益 曼 法 模拟 的 植保 无 人 机 施 药 地 面 沉 


梭 与 田间 实验 结果 对 比 '” 


0 2 
展 向 距离 X /m 


Fig. 6 Droplet deposition distributions of UAV spraying 


based on LBM simulation and field test 


积分 布 等 参数 的 实际 测量 予以 验证 。 目 前 主要 验 
证 手段 主要 有 田间 采样 和 室内 测量 两 种 。 


3.1 田间 采样 方法 


传统 航空 施 药 实 验 主要 采用 田间 试验 手段 ， 
一 般 利 用 离散 式 地 面 雾 滴 沉 积 样本 检测 及 空间 离 
散 式 风 速 测量 数据 为 依据 对 数值 模拟 结果 进行 评 


ChinaXiv 合 作 期 刊 


202302.00201v1 


chinaXiv 


fh ON ， 并 采用 高 速 风 洞 试验 获得 初始 喷雾 粒 
径 分 布 信息 “*”。 而 植保 无 人 机 旋 流 下 洗 气 流 复 
杂 ， 田 间 离 散 式 风 速 测量 难以 获得 瞬时 流动 结 
构 ， 且 其 距离 靶 标 近 ， 喷 雾 初 始 速度 对 实际 筋 滴 
沉积 的 影响 较 大 。 传 统 测量 手段 在 处 理 植 保 无 人 
机 施 药 过 程 测量 方面 面临 一 定 困难 。 

在 速度 场 测量 方面 ， 主 要 通过 点 阵 分 布 式 测 
量 手 段 大 规模 布设 速度 测 点 ， 进 行 植保 无 人 机 旋 
辟 下 洗 风 场 测量 。 如 李 继 宇 等 ”利用 田间 无 线 
微 气象 风速 传 感 带 阵列 测量 多 旋 注 植保 无 人 机 下 
洗 风 场 法 ， 发 现实 际 作 业 风 场 存 在 风 幕 效应 ， 风 
幕 两 侧 风 速 存 在 突变 。 王 昌 陵 等 ” 同样 基于 田 
间 无 线 微 气象 风速 传感器 阵列 获取 了 三 种 不 同 植 
保 无 人 机 在 不 同方 向 上 的 下 洗 气 流速 度 场 分 布 。 
杨 风 波 等 ”采用 转 轮 式 风速 计 逐 点 测量 多 旋 恤 
植保 无 人 机 悬 停 状 态 下 的 风 场 分 布 情况 ， 用 于 验 
证 数值 模型 对 下 洗 风速 的 模拟 精度 ， 实 验 测量 结 
果 与 模拟 结果 差异 在 10% 以内。 上述 研 究 采用 的 
转 轮 式 传感器 成 本 较 低 ， 适 合 进行 阵列 式 布置 开 
展 大 范围 田间 实验 ， 但 由 于 存在 转子 旋转 惯性 ， 
其 对 高 频 气 流 方向 变化 的 响应 频率 较 低 ， 在 测量 
植保 无 人 机 旋 避 下 洗 气 流 这 种 非 定常 流 场 时 效果 
较 差 ， 在 气流 速度 变化 幅度 明显 的 区 域 延迟 时 间 


ma f d ca 
* fab = 


(a) ELE XL d 


较 长 ， 影 响 测 量 结果 准确 性 。 为 克服 这 一 问题 ， 
其 他 更 高 频 测 量 手段 相继 被 应 用 。 如 Zhang 等 
采用 热线 风速 仪 逐 点 测量 多 旋 副 植保 无 人 机 其 停 
状态 下 的 风 场 分 布 情况 ， 用 于 验证 其 数值 模型 模 
拟 准 确 性 ， 实 测 结果 与 数值 模拟 结果 误差 低 于 
10%。 热 线 风 速 仪 具备 1 kHz 以 上 的 测量 频率 ， 
适合 于 对 非 定常 流 场 开展 瞬时 高 精度 测量 ， 但 由 
于 成 本 较 高 ， 一 般 不 采用 阵列 式 布置 ， 仪 能 通过 
单 点 采样 获得 整体 风 场 中 的 速度 及 脉动 量 平均 分 
布 。 因 此 该 方案 目前 还 主要 用 于 固定 高 度 的 植保 
无 人 机 及 停 风 场 测量 中 ,保持 各 项 条 件 不 变 的 状 
态 下 改变 测 点 位 置 ， 以 提供 空间 平均 场 测 量 
结果 。 

针对 上 述 方法 各 目的 局 限 性 ， 作 者 团队 采用 
三 维 超声 风速 仪 阵 列 开展 航空 施 药 风 场 测量 人 研 
究 "。 三 维 超声 风速 仪 测量 频率 达 20 Hz， 进 行 
空间 阵列 式 布置 后 ， 能 够 捕捉 部 分 高 频 涡 变 化 。 
作者 团队 分 别 测量 了 有 人 驾驶 固定 翼 飞 机 尾 涡 运 
动 耗 散 过 程 ， 以 及 植保 无 人 机 悬 停 及 巡航 作业 状 
d PAY PUE ULM, ARBEIT E FUN 
内 测量 及 数值 模拟 奠定 了 基础 ， 测 试 场景 如 图 7 
所 示 。 


(b) 植 保 无 人 机 
图 7 运用 三 维 超声 风速 仪 阵列 测量 固定 翼 飞 机 和 植保 无 人 机 瞬时 风 场 


Fig. 7 Instantaneous flow field measurement of fixed-wing aerial plane and UAV with 3D ultrasonic anemometer array 


雾 滴 沉 积分 布 测量 方面 ， 传 统 水 敏 纸 或 聚 乙 
类 卡片 等 单 点 式 地 面 采 样 方式 依然 是 目前 植保 无 
人 机 施 药 药 液 沉 积 允 移 实验 最 主流 的 测量 方 
式 “。 但 由 于 植保 无 人 机 下 洗 风 场 的 非 定常 特 
性 ， 其 造成 的 雾 滴 沉积 分 布 呈 现 非 均匀 特性 ， 单 


点 式 采 样 手段 受 空间 分 辩 率 限制 ， 准 确 性 较 差 。 
作者 团队 °°" 针对 上 述 问 题 ， 开 发 了 连续 式 荧光 
纸 带 采样 仪 ， 可 在 喷洒 幅 宽 内 进行 连续 式 采 样 ， 
从 而 获得 高 空间 分 辩 率 的 雾 滴 沉 积 量 分 布 结果 。 
其 测量 及 分 析 过 程 如 图 8 所 示 。 
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图 8 基于 光谱 分 析 的 航空 施 药 雾 滴 沉 积 特性 检测 系统 


Fig. 8 Aerial spray droplet deposition test system based on spectrum analysis 


男 外 ， 大 范围 非 接触 式 测量 如 红外 热 成 像 
法 、 激 光 雷 达 、 近 红外 开路 傅 里 叶 变换 方法 等 相 
继 出 现 ， 一 定 程 度 上 增加 了 筋 滴 运 动 田间 实验 方 
式 的 多 样 性 ， 提 升 了 雾 滴 运 动 过程 的 直观 性 理 
解 。 例 如 ， 张 京 等 UU 采用 红外 热 像 仪 与 无 人 机 
联 用 ， 研 究 了 无 人 驾驶 直升机 航空 喷雾 参数 对 药 
液 筋 滴 沉 积 效果 的 影响 。 通 过 测试 喷雾 前 、 后 作 
物 冠 层 温 度 变化 率 反 映 筋 滴 在 水 稻 冠 层 的 沉积 效 
果 ， 确 定 红 外 热 成 像 技 术 可 以 准确 反映 雾 滴 在 水 
稳 上 的 沉积 规律 。Jiao 等 ”也 运用 红外 热 像 仪 
对 固定 辟 飞 机 施 药 筋 滴 云 团 在 空间 中 的 运动 沉积 
过 程 进 行 了 记录 ， 成 功 提取 了 雾 滴 云 团 的 沉降 过 
程 及 影响 范围 。Gregorio $ "' 采用 激光 雷达 扫 
描 喷 雾 团 ， 获 取 具 有 时 间 和 空间 分 辨 率 的 数据 。 
激光 雷达 喷雾 味 移 测量 值 与 水 平 采样 点 获得 的 测 
量 值 比 较 ， 其 决定 系数 高 达 0.85 以 上 。 激 光 雷 达 
系统 还 可 用 于 监测 喷雾 雾 团 的 演化 过 程 并 生成 这 
些 雾 团 的 二 维 图 像 。Kira 等“ 利用 近 红 外 开路 
傅 里 叶 变 换 方法 监测 及 表征 药 液 筋 滴 肚 移 问题 ， 
研究 发 现在 采用 水 敏 纸 估计 药 液 雾 滴 粒 径 分 布 
后 ， 该 方法 能 够 在 扫描 线 内 获得 筋 滴 载 量 ， 同 时 
光谱 特征 还 能 用 于 鉴定 农药 溶液 的 有 机 相 。 

但 上 述 方法 仍 存在 成 本 较 高 、 空 间 分 辩 率 
低 、 受 环境 复杂 背景 干扰 、 空 间 积 分 效应 难以 去 


除 等 问题 ， 相 对 传统 采样 方法 ， 仍 未 能 大 规模 应 
用 于 航空 施 药 田间 实验 中 。 


3.2 室内 测量 方法 


由 于 田间 实验 受 环 境 影 响 因 素 较 大 ， 日 实际 
作业 区 域 大 ， 部 分 精细 化 测量 手段 难以 实现 ， 因 
此 ， 室 内 测量 方法 作为 田间 实验 的 补充 ,适用 于 
获得 更 多 非 定常 流 场 的 验证 数据 ， 可 验证 不 同 数 
值 模拟 方法 对 复杂 旋 波 风 场 模拟 效果 。 目 前 可 用 
于 室内 测量 植保 无 人 机 下 洗 风 场 和 雾 滴 运 动 的 方 
法 主要 有 热线 风速 仪 、 激 光 雷 达 、 激 光 粒 度 仪 、 
PIV、 相 位 多 普 勒 测速 (Phase Doppler Interfer- 
ometry, PDI) 等 。 由 于 田间 实验 测量 方法 中 已 
介绍 了 热线 风速 仪 、 激 光 雷 达 等 测量 方法 ， 本 节 
拟 重点 介绍 激光 粒度 仪 、PIV 和 PDI 系 统 。 

3.2.1 激光 粒度 仪 

激光 粒度 仪 基于 激光 衍射 原理 (Laser Dif- 
fraction), ， 可 对 各 种 物料 和 样品 进行 粒度 检测 , 
并 在 较 短 的 时 间 内 给 出 比较 详细 的 粒度 分 布 数 
据 ， 适 用 于 在 室内 开展 喷雾 粒 径 谱 测 量 工 作 。 

该 仪 需 利用 光 的 衍射 现象 ， 即 不 同 粒 径 颗 粒 
产生 的 衍射 角 差 异 ， 通 过 计算 探测 器 上 收集 到 的 
不 同 衍射 图 形 的 光 强 分 布 ， 给 出 颗粒 的 粒度 大 小 
和 粒度 分 布 。 不 同 大 小 的 颗粒 的 衍射 光 强 集中 在 
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探测 右 的 不 同 部 位 ， 根 据 在 多 元 探测 带 上 得 到 的 
衍射 光 强 的 分 布 ， 通 过 颗粒 大 小 和 光 强 分 布 之 间 
的 相关 公式 计算 得 到 颗粒 的 粒度 分 布 。 其 典型 仪 
器 系统 及 测量 原理 见 图 9。 


接收 器 â 


(a) 激 光 粒 度 仪 典型 样机 


ze 1555; "E 1-3 i ea 
BOLE HOE LES ORE 接收 器 


wae 


(pb) 激 光 粒 度 仪 喷雾 粒 径 测量 原理 
图 9 激光 粒度 仪 及 其 喷雾 粒 径 测量 原理 示意 
Fig.9 Schematic diagram of the laser diffraction system and 


the droplet diameter measurement 


3.2.2 PIV 系统 

PIV 系统 包括 激光 器 、 片 光源 成 型 器 、 相 
机 、 同 步 器 ， 及 粒子 图 像 测速 软件 等 ， 基 于 激光 
平面 粒子 成 像 ， 通 过 对 粒子 图 像 时 间 序 列 进行 相 
关 计 算 ， 从 而 获得 流动 瞬时 速度 场 ， 其 典型 仪器 
系统 如 图 10 所 示 。 


计算 软件 


同步 控制 器 
图 10 典型 粒子 图 像 测速 系统 


Fig. 10 Typical particle image velocimetry system 

采用 高 速 激光 器 及 高 速 互 补 金属 氧化 物 半 导 
体 (Complementary Metal Oxide Semiconductor, 
CMOS) 相机 的 时 间 解 析 PIV 技术 ， 是 目前 唯一 
能 够 获得 无 人 机 下 洗 风 场 瞬时 流动 结构 的 机 理性 


目前 该 仪 吉 已 成 为 各 喷雾 测量 实验 室 的 标准 
配置 ， 如 美国 农业 部 农业 服务 署 (USDA-ARS), 
澳大利亚 昆士兰 大 学 (The University of 
Queensland) 、 美 国内 布 拉 斯 加 大 学 (The Univer- 
sity of Nebraska) 、 中 国 农业 大 学 、 北 京 农业 智能 
装备 技术 研究 中 心 等 单位 均 配 备 了 该 仪器 ， 主 要 
针对 地 面 /航空 施 药 用 喷头 开展 喷雾 粒 径 分 布 测 
量 实验 ， 获 得 不 同类 型 喷头 喷雾 粒 径 分 布 结 
果 ““ ， 基 于 该 仪器 获取 并 积累 的 大 量 数据 也 促 
成 了 地 面 /航空 施 药 喷头 分 类 ASABE 及 ISO 标准 
的 建立 。 

激光 粒度 仪 雾 滴 粒 径 测量 范围 广 ， 可 达 
0.1—3500 pgm， 扫描 频率 可 达 2000 Hz， 适 用 于 
对 高 密度 雾 团 进行 快速 扫描 检测 。 但 其 功能 
一 ， 只 能 测量 雾 滴 粒 径 谱 ， 且 只 能 对 激光 路 径 上 
的 全 部 雾 滴 积分 测量 ， 无 法 获得 雾 滴 粒 径 三 维 空 
间 分 布 的 问题 。 


研究 手段 。 在 无 人 机 下 洗 风 场 测量 中 ， 能 够 提供 
mA RMA, SCL ie EBL 
AY EFE Wet, CDE TR RT ER Ts AD 
准确 度 起 决定 性 作用 ， 如 图 11 所 示 。 

目前 该 领域 国际 领先 的 研究 机 构 包 括 德国 德 
宇航 中 心 (Deutsches Zentrum für Luftund 
Raumfahrt, DLR), 3) 5 H ŽK% (Universi- 
ty of Maryland) 以 及 新 加 坡 国 立 大 学 (National 
University of Singapore) 等 。 其 中 德国 DLR 的 
Raffel 是 PIV 技术 发 展 的 主要 贡献 者 ， 领 导 了 欧 
洲 直 升 机 实验 测试 项 目 ， 其 团队 运用 PIV 技术 对 
旋 疲 己 面 、 机 身 附 近 流 场 进行 了 细致 测量 ， 主 要 
用 于 分 析 不 同 流 态 下 旋翼 翼 面 效率 和 受 力 等 问 
题 ， 较 少 涉 及 旋翼 下 洗 气 流 和 周边 气流 的 研 
FTO I E RH Leishman IBA CE fe Se XU 
测量 方面 具有 多 年 研究 经 验 ， 综 合 运 用 数值 模 
拟 、 流 动 显示 技术 、 高 速 PIV 测量 技术 等 ， 对 单 
旋翼 下 洗 风 场 的 精细 结构 发 展 规律 进行 了 深入 研 
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(b) 基 于 LBM 模 拟 的 旋翼 翼 炎 涡 涡 量 分 布 


图 11 基于 PIV 试验 结果 验证 和 LBM 对 姻 尖 涡 耗 散 过 程 模拟 的 准确 性 


Fig. 11 Accuracy verification with PIV experiment to the and LBM modeling of rotor tip vortex decay process 


究 。 但 其 局 限于 室内 实验 条 件 ， 将 研究 重点 放 在 
微型 旋 副 研究 上 ， 其 研究 方法 可 以 借鉴 ， 但 研究 
成 果 难 以 直接 复制 到 植保 无 人 机 作业 中 '“”。 新 
加 坡 国立 大 学 的 Nathan 则 运用 动 壁面 和 高 速 PIV 
技术 研究 了 旋 恤 在 不 同 前 飞速 度 下 的 风 场 演化 过 
程 ， 该 实验 考虑 了 前 飞 状 态 下 旋翼 与 地 面 的 相对 
运动 对 下 洗 气 流 结构 演化 的 影响 ， 但 仍然 局 限于 
室内 实验 的 模型 尺度 限制 UU. EXRÜP LAUDE 
无 人 机 非 定 常 旋 转 流 场 的 测量 方面 开展 了 大 量 基 
础 研究 工作 ， 但 受 室 内 实验 条 件 限 制 ， 尚 未 涉及 
包含 无 人 机 机 体 的 大 范围 流 场 测量 研究 。 针 对 上 
述 情况 ， 作 者 团队 率先 针对 AF25B AICHE MGIP 
等 植保 无 人 机 开展 相关 室内 缩 比 模型 PIV 瞬时 流 
yw HF, AAG TAKER RA BRE BGT 
程 、 典 型 涡 结构 、 平 均 速 度 场 及 流 线 图 等 结果 ， 
并 与 团队 先期 进行 的 田间 作业 风 场 测试 结果 和 开 
发 的 该 型 无 人 机 数值 模拟 结果 分 别 进行 对 照 ， 构 
建 了 从 田间 实验 到 数值 模拟 的 联接 桥梁 。 
3.23 PDI 系统 

PIV 系统 主要 针对 植保 无 人 机 旋 辟 下 洗 风 场 
发 展演 化 过 程 进行 研究 ， 而 对 植保 无 人 机 施 药 药 
液 初始 运动 过 程 则 主要 依靠 PDI 系 统 进 行 研究 。 
由 于 植保 无 人 机 施 药 过 程 距 离 冠 层 近 ， 初 始 喷雾 
雾 滴 速度 分 布 对 沉积 分 布 的 影响 不 可 忽略 ， 传 统 


用 于 雾 滴 粒 径 测 量 的 激光 粒度 仪 等 装备 无 法 满足 
需求 。PDI 技 术 运 用 二 极 管 泵 浦 固体 激光 器 构成 
的 光学 系统 ， 采 用 激光 相位 多 普 勒 效应 ， 可 同时 
测量 通过 激光 焦点 的 雾 滴 粒 径 及 速度 。 其 粒 径 测 
量 范围 为 0.3~7000 hum， 速度 测量 范围 为 -100~- 
300 m/s， 激 光 相 位 多 普 勒 测速 仪 见 图 12。 


n es 
图 12 激 位 多 首 勒 测速 系统 
Fig.12 Phase Doppler interferometry system 

该 方法 可 应 用 于 各 型 喷头 的 喷雾 粒 径 及 速度 
分 布 测量 ， 如 Nuyttens 4& 7577 对 不 同类 型 喷头 喷 
筋 区 的 粒 径 及 速度 分 布 进行 了 细致 测量 ， 分 析 了 
其 对 雾 滴 沉积 球 移 的 潜在 影响 。Li 等 7" 针对 低 
压条 件 下 的 平面 扇形 喷头 ， 对 在 不 同 高 度 下 喷雾 
书面 椭圆 截面 的 粒 径 及 速度 分 布 进行 了 统计 ， 基 
于 喷雾 角度 、 流 量 、 孔 口 直径 等 参数 ， 获 得 了 雾 
滴 粒 径 拟 合 方程 。 已 有 研究 尚未 开展 复杂 风 场 作 
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用 下 的 雾 滴 粒 径 及 速度 分 布 变 化 规律 ， 其 初始 速 
度 分 布 测量 结果 也 未 作为 参考 条 件 被 应 用 于 数值 
模拟 中 。 

作者 团队 “和 针对 上 述 研究 现状 ， 基 于 PDI 
测试 平台 测量 了 植保 无 人 机 常用 平面 扇形 噶 头 产 
生 的 初始 雾 滴 粒 径 及 速度 分 布 ， 结 合 数 值 模型 分 
析 了 其 在 旋 恤 风 场 作用 下 的 粒 径 及 速度 空间 分 布 
变化 情况 ， 为 实现 雾 滴 运 动 沉 积 球 移 预测 提供 了 
初始 喷雾 参数 支撑 。 


4 综合 分 析 


目前 无 粘 模型 、 计 算 流 体力 学 模型 《有限 体 
积 法 及 有 限 差 分 法 ) 及 格子 玻 尔 效 曼 模型 在 模拟 
植保 无 人 机 施 药 过 程 中 均 有 应 用 ， 且 具备 各 上 自 的 
优势 、 缺 陷 及 适用 范围 ， 各 模型 适用 性 、 复 杂 度 
等 优 缺点 比较 见 列 表 2。 

表 2 各 数值 模型 针对 植保 无 人 机 施 药 模 拟 的 优 缺 点 比较 


Table 2 Comparison of advantages and disadvantages of nu- 


merical models for plant protection UAV spraying 


冠 层 适 环境 适 模型 适 消 流 适 计算 量 操作 涡 模 拟 
y 


模型 ”用 性 ”用 性 ”用 性 ”用 性 难度 ”精度 
无 粘 模型 低 中 中 ik tk od 中 
有 限 体积 法 中 高 中 高 高 中 中 
有 限 差分 法 中 高 低 高 极 高 高 高 


A 高 低 高 中 高 低 sm 


无 粘 模 型 模拟 植保 无 人 机 施 药 过 程 拥 有 计算 
资源 消耗 低 ， 对 近 场 范围 沪 尖 涡 发 展 模拟 准确 度 
较 好 ， 远 场 有 人 机 施 药 雾 滴 肚 移 模拟 准确 性 经 过 
大 量 验证 ， 适 用 复杂 条 件 等 优势 。 其 面临 的 主要 
问题 是 涡 元 法 中 疲 尖 涡 结构 不 耗 散 的 特性 在 远 场 
不 符合 实际 情况 ， 目 前 暂时 无 法 处 理 大 空间 范围 
的 模拟 问题 。 男 外 涡 元 法 采用 的 三 维 模 拟 和 二 维 
雾 滴 远 场 肚 移 模型 的 结合 目前 仍 存在 衔接 难点 ， 
尚 无 相关 文献 验证 其 结合 效果 。 二 维 雾 滴 味 移 模 
型 中 的 简化 冠 层 模型 也 无 法 用 于 涡 元 法 三 维 模型 
中 ， 进 一 步 制 约 了 无 粘 模 型 在 模拟 植保 无 人 机 在 
包含 冠 层 影响 的 实际 作业 条 件 下 的 模拟 能 

计算 流体 力学 模型 中 ， 有 限 体积 法 基础 完 
备 ， 适 应 范围 广 ， 可 扩展 性 强 ， 是 主流 商用 软件 


的 首选 方法 ， 适 合 对 复杂 和 气象、 地 形 条 件 下 的 植 
保 无 人 机 作业 引起 的 风 场 变化 和 雾 滴 远 场 际 移 影 
响 进 行 模拟 。 其 主要 面临 格式 精度 低 易 造成 人 工 
耗 散 、 非 结构 网 格 条 件 下 的 网 格 自 适应 质量 难以 
保证 、 大 范围 三 维 模 拟 计算 量 过 大 等 问题 。 由 于 
有 限 体积 法 格式 精度 较 低 ， 如 果 采 用 正常 计算 量 
且 不 加 入 自 适 应 网 格 , 会 出 现 波 尖 涡 耗 散 速 度 远 
快 于 实际 情况 的 问题 ， 从 而 影响 风 场 及 雾 滴 运 动 
模拟 准确 度 。 在 目前 与 植保 无 人 机 相关 的 模拟 研 
究 中 ， 该 现象 普遍 存在 。 而 要 准确 模拟 辟 尖 涡 运 
动 耗 散 过 程 ， 需 运用 更 多 的 自 适应 网 格 来 细 化 副 
尖 涡 发 展区 域 ， 从 而 急剧 增加 模拟 计算 量 。 在 模 
拟 大 范围 雾 滴 丈 移 过 程 时 ， 计 算 量 消耗 问题 依然 
是 该 方法 的 瓶颈 。 有 限 差 分 法 可 轻易 构造 高 精 格 
式 ， 适 用 于 简单 外 形 模型 产生 的 精细 风 场 结构 模 
拟 ， 目 前 主要 应 用 于 基础 研究 领域 ， 工 程 应 用 较 
少 。 其 在 解决 复杂 外 形 的 结构 化 网 格 构建 方面 依 
然 存在 困难 ， 尤 其 是 面 对 无 人 机 复杂 结构 外 形 及 
冠 层 复杂 结构 时 。 这 进一步 限制 了 该 方法 在 植保 
无 人 机 风 场 与 冠 层 作用 模拟 方面 的 应 用 。 

基于 LBM 的 数值 模型 具备 四 阶 模 拟 精度 、 
自 适应 网 格 及 网 格 自动 生成 功能 ， 适合 于 快速 对 
复杂 外 形 植保 无 人 机 及 冠 层 作用 下 的 风 场 分 布 及 
雾 滴 运动 进行 模拟 。 其 主要 问题 是 固有 立方 体 网 
格 在 需要 进行 复杂 外 形 边界 层 加 密 时 网 格 数 急剧 
上 升 。 为 兼顾 计算 量 适当 降低 边界 层 网 格 数 ， 则 
会 影响 模型 模拟 精度 ， 尤 其 会 对 植保 无 人 机 旋 波 
造成 的 下 洗 速 度 值 产生 较 大 影响 。 男 外 该 模型 在 
雾 滴 运动 模拟 方面 仍 基 于 初级 的 受 力 模型 ， 尚 未 
涉及 环境 温 湿度 对 雾 滴 莹 发 的 影响 ， 目 前 仍 无 法 
模拟 复杂 环境 条 件 下 的 植保 无 人 机 施 药 过 程 。 

在 数值 模型 验证 方法 方面 ， 目 前 田间 实验 
(超声 / 风 杯 风速 仪 ， 水 敏 纸 及 聚 乙烯 卡片 ， 荧 光 
纸 带 ， 非 接触 测量 方法 ) 及 室内 实验 (PIV, 
PDI) 等 验证 方法 在 验证 植保 无 人 机 施 药 数值 模 
型 模拟 准确 性 过 程 中 均 有 应 用 ， 且 具备 各 自 的 优 
势 、 缺 陷 及 适用 范围 ， 各 方法 优 缺 点 比较 见 
列表 3。 
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表 3 各 数值 模型 验证 方法 的 优 缺 点 比较 
Table 3 Comparison of advantages and disadvantages of verification methods for numerical model 

模型 验证 方法 风速 雾 滴 速 度 /位 置 ” mE 时 间 分 辨 灰 ”空间 分 辨 率 ”场地 适应 性 操作 难度 
超声 风速 仪 Y -=> -= 高 低 高 中 
风 杯 式 风 速 仪 Y — AER 低 低 低 
水 敏 纸 及 卡片 = — Al su 中 高 低 
荧光 纸 带 I - J "me = a E 
热线 风速 仪 Y —— —— 高 低 中 高 
红外 热 成 像 EE y = 中 中 高 中 
激光 雷达 一 一 高 d 高 " 
近 红 外 开路 傅 里 叶 Y - -一 中 中 高 中 
激光 粒度 仪 == —— Y 中 中 中 中 
PIV y Y y 高 高 低 高 
PDI 一 一 Y y 高 中 低 高 


超声 风速 仪 及 风 杯 式 风速 仪 是 测量 植保 无 人 
机 田间 作业 风 场 最 重要 的 手段 ， 可 用 于 验证 数值 
模拟 风 场 速度 准确 性 。 超 声 风速 仪 具有 高 测量 频 
率 ， 但 价格 昂贵 ， 体 积 较 大 且 为 单 点 式 测 量 ， 使 
用 普及 性 较 低 。 风 杯 式 风速 仪 体 积 较 小 ， 成 本 
1K, 适合 大 范围 布置 但 测量 响应 频率 很 低 ， 无 
法 处 理 风 场 波动 幅度 较 大 的 情况 。 

水 敏 纸 、 肾 乙烯 卡片 及 区 光 纸 带 主要 用 于 测 
量 植保 无 人 机 施 药 幅 宽 内 药 液 沉积 分 布 情况 ， 可 
用 于 验证 数值 模型 对 雾 滴 地 面 沉积 分 布 预测 准确 
性 。 其 中 水 敏 纸 及 聚 乙 烦 卡片 采用 离散 式 采样 布 
置 方 式 ， 其 采样 空间 分 辨 率 受 限于 布点 密度 ， 尤 
其 是 对 植保 无 人 机 作业 风 场 这 一 非 定常 流 场 导致 
的 空间 分 布 波动 极 大 的 雾 滴 沉 积 情况 难以 真实 反 
WR DEI AEC A BP IE BER EE Ta UK Hb T ERA 
E, RAKEI TSF Te TUR AR FE Ze le] PRS, ROS 
较 好 的 应 用 前 景 。 

大 范围 非 接触 式 测量 如 红外 热 成 像 法 ， 激 光 
雷达 ， 近 红外 开路 傅 里 叶 变 换 方 法 等 能 够 在 不 干 
扰 被 测 对 象 的 前 提 下 提供 具备 时 间 、 空 间 分 辨 率 
的 测量 结果 ， 提 升 了 航空 施 药 雾 滴 运 动 过 程 的 直 
观 性 理解 。 但 上 述 方法 仍 存在 成 本 较 高 、 空 间 分 
辨 率 低 、 受 环境 复杂 背景 干扰 、 空 间 积 分 效应 难 
以 去 除 等 问题 ， 相 对 传统 采样 方法 ， 仍 未 能 大 规 
模 应 用 于 航空 施 药 田间 实验 中 。 

激光 粒度 仪 是 各 喷雾 测量 实验 室 的 标准 配 
置 ， 具 有 测量 简便 快速 的 优势 。 但 其 功能 单一 ， 


只 能 测量 激光 路 径 上 所 有 筋 滴 粒 径 的 平均 分 布 ， 
无 法 完全 满足 室内 实验 的 多 样 化 需求 。 而 室内 
PIVPDI 实 验 主要 作为 田间 实验 的 补充 ， 可 对 难 
以 开展 田间 测量 的 植保 无 人 机 瞬时 下 洗 风 场 ， 喷 
雾 粒 径 及 速度 分 布 等 信息 进行 测量 ， 和 田间 实验 
结果 一 起 共同 为 植保 无 人 机 施 药 数值 模型 提供 验 
证 信息 。 其 中 PIV 因 其 具备 高 时 空 分 辨 率 ， 以 及 
IET ER BRE HE TE, EEI AKH A IR EHE 
粒 径 等 信息 进行 测量 ， 但 其 对 操作 要 求 较 高 ， 环 
境 适 应 性 不 强 。PDI 可 对 单 点 进行 速度 和 粒 径 扫 
描 测 量 ， 其 面临 的 主要 问题 是 时 空 分辨 率 不 能 和 
顾 ， 一 般 用 于 对 喷头 初始 喷雾 特性 进行 测量 ， 为 
数值 模型 提供 喷雾 初始 条 件 。 


5 展望 


随 着 研究 的 不 断 推进 ， 关 于 植保 无 人 机 施 药 
雾 滴 沉 积 味 移 研 究 的 各 类 模拟 方法 未 来 将 继续 改 
进 优化 。 上 述 三 种 模拟 方法 将 根据 自身 特性 ， 适 
应 不 同 的 研究 需求 。 例 如 一 般 工程 应 用 中 要 求 快 
速 运算 ， 因 此 仍 以 无 粘 横 型 为 主 ; 涉及 流动 细节 
分 析 和 喷 施 优化 的 研究 则 越 来 越 多 采用 格子 玻 尔 
效 曼 法 ; 对 于 复杂 气象 条 件 、 基 础 性 流动 结构 及 
两 相 流 的 研究 ， 计 算 流 体力 学 模型 依然 是 最 佳 
选择 。 

无 粘 模型 中 ， 近 场 三 维 模 拟 与 二 维 远 场 雾 滴 
味 移 模拟 两 种 方法 间 的 接口 将 继续 优化 ， 最 终 形 
成 统一 的 植保 无 人 机 快速 模拟 方案 ， 同 时 进一步 
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发 挥 该 方法 的 计算 速度 优势 ， 向 基于 植保 无 人 机 
飞行 实测 数据 的 施 药 雾 滴 沉 积 球 移 实时 模拟 演 
进 。 该 模型 可 能 会 成 为 植保 无 人 机 作业 地 面 站 的 
标准 配置 ， 能 够 根据 作业 及 环境 参数 实时 预测 植 
保 无 人 机 施 药 作 业 过 程 中 的 雾 滴 潜 在 球 移 区 ， 从 
而 有 针对 性 地 重新 规划 路 线 或 变更 设置 作业 
条 件 。 

计算 流体 力学 模型 中 ， 有 限 体积 法 在 自 适应 
网 格 及 并 行 算法 的 支持 下 ， 其 人 工 耗 散 率 高 ， 计 
算 开销 大 的 问题 能 获得 一 定 的 缓解 。 在 今后 依然 
是 模拟 复杂 气象 、 作 业 条 件 下 植保 无 人 机 风 场 及 
雾 滴 运 动 过 程 并 获取 精确 结果 的 最 佳 手段 。 该 方 
法 适合 于 针对 植保 无 人 机 机 型 38 A S 4 ces 
调整 带 来 的 气动 特性 和 风 场 分 布 变化 进行 机 理 
HE, AEW. KEZMA., Rmi, K 
多 种 自 定义 函数 功能 ， 本 方法 还 适合 对 复杂 气 
象 、 地 形 条 件 下 的 植保 无 人 机 作业 引起 的 风 场 变 
化 和 雾 滴 远 场 球 移 影 响 进 行 模拟 。 有 限 差 分 法 未 
来 将 继续 发 挥 其 模拟 精度 优势 ， 在 基础 研究 领域 
对 特定 机 型 产生 的 尾 涡 结构 进行 精细 化 模拟 ， 与 
室内 实验 结果 互 为 补充 ， 为 其 他 数值 模型 提供 更 
为 精确 的 尾 涡 运动 耗 散 过 程 对 照 ， 帮 助 验证 和 提 
升 其 他 模型 模拟 准确 度 。 

基于 LBM 的 数值 模型 在 对 植保 无 人 机 作业 
风 场 模拟 、 雾 滴 运 动 及 球 移 过 程 模 拟 施 药 参 数 、 
药械 结构 、 喷 头 选择 优化 等 方面 都 具有 不 可 替代 
的 作用 ， 适 合 于 快速 对 复杂 外 形 植保 无 人 机 及 冠 
层 作 用 下 的 风 场 分 布 及 雾 滴 运 动 进行 模拟 。 该 模 
型 未 来 需要 进一步 优化 其 壁 滑 流 模型 以 克服 粗 网 
格 在 壁面 附近 造成 的 模拟 精度 损失 ， 平 衡 计算 量 
与 计算 精度 的 关系 。 同 时 需要 进一步 加 强 雾 滴 蒸 
发 、 碰 撞 、 破 碎 模 型 及 环境 温 湿度 模型 。 

数值 模型 验证 方法 方面 ， 田 间 实 验 将 继续 保 
持 特性 ， 并 弥补 受 环境 影响 较 大 、 时 空 分 辨 率 不 
高 的 问题 。 主 要 包括 : (1) 风 场 测量 环节 发 展 廉 


价 、 小 型 化 三 维 超声 风速 仪 及 热线 风速 仪 等 高 时 
间 分 辩 率 测量 方法 ， 并 缩小 测 点 间隔 ， 增 加 测 点 


数量 以 弥补 空间 分 辨 率 不 足 ; (2) 雾 滴 沉 积 检测 


环节 发 展 自动 化 荧光 纸 带 检 测 等 连续 性 采样 方 
法 ,提升 空 间 分 辨 率 ; (3) 环境 背景 布置 高 频 风 
场 测量 设备 ， 同 步 扣 除 背 景 噪 音 以 抵消 环境 风 场 
变化 对 测量 结果 造成 的 随机 性 影响 ; (4) 雾 滴 运 
动 /位 置 测定 主要 发 展 非 接 触 式 测量 手段 ， 如 红 
外 热 成 像 、 激 光 雷 达 、 近 红外 开路 傅 里 叶 方 法 
等 ， 着 力 提 高 其 时 空 分辩 率 及 抗 环境 背景 干扰 
能 力 。 

室内 实验 除 传统 激光 粒度 仪 测量 方式 外 ， 将 
以 PIV 和 PDI 等 高 频 测 量 手段 为 基础 ， 一 方面 继 
续 提 升 系统 自动 化 及 智能 化 水 平 ， 降 低 操作 难度 
和 布置 要 求 ， 提 升 室 内 测量 效率 。 一 方面 继续 拓 
展 应 用 的 适用 范围 ， 例 如 积极 拓展 室外 大 范围 
PIV SG & 25 ig WE LN. fes PDI My Se SH 
测量 系统 等 ， 实 现 室内 外 测量 结果 的 自然 衔接 , 
为 数值 建 模 提 供 更 全 面 的 验证 信息 支持 。 
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Abstract: With the increasing application of plant protection unmanned aerial vehicle (UAV) in precision agriculture, the numeri- 
cal simulation methods for the development of the downwash flow field of the plant protection UAV and the deposition and drift 
process of droplets affected by the downwash flow field have achieved rapid and diversified development, but the advantages, 
disadvantages, scope of application, and verification of each method still lack a systematic review. This article discusses the in- 
viscid model, computational fluid dynamics model and lattice Boltzmann model (LBM) respectively. The advantage of the invis- 
cid wake vortex model based on the vortex element method is that the calculation process is simple. Moreover, integrated with 
the most widely used aerial spray drift prediction software AGricultural DISPersal (AGDISP), it can be a promising way to do 
real-time UAV spray drift prediction. But due to lack of viscosity and turbulence models, the droplet deposition and drift simula- 
tion accuracy of inviscid model is relatively lower than other models. The computational fluid dynamics (CFD) model includes 
the finite volume method (FVM) and the finite difference method (FDM). The FVM in the computational fluid dynamics model 
has high robustness and can be applied to the simulation of various complex environments. Many commercial CFD software are 
based on FVM and achieved a fast development in aerial spray modeling recently. However, the FVM is greatly affected by the 
quality of the mesh, and its commonly used upwind style has limited accuracy (second-order accuracy). Under the same mesh 
density, it is easier to generate artificial dissipation when simulating the rotor tip vortex than the finite difference method. As a 
result, the simulated rotor tip vortex dissipation speed is much faster than the actual situation. Compared with the FVM, the 
structured grid used in the FDM is easier to construct a high-order precision numerical format. Which can reach 4-5 orders of 
accuracy, and with adaptive grid technology, FDM can simulate the evolution of rotor tip vortex with high temporal and spatial 
accuracy, and can reproduce the typical flow structure development process of the real rotor downwash flow field. However, it 
also has problems such as high grid structure requirements and excessive computing power requirements. LBM has advantages 
in computing three-dimensional flow field problems with complex boundary conditions and non-stationary moving objects. 
However, there are still shortcomings in its functional diversity and completeness. The accuracy of the numerical models men- 
tioned above still needs field test and indoor experiment such as high-speed Particle Image Velocimetry (PIV)/ Phase Doppler 
Interferometry (PDI) method to verify and optimize. The authors finally pointed out the future direction of plant protection UAV 
application simulation and verification. 
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